Ion-acoustic solitary waves and shocks in a collisional dusty negative ion plasma 
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We study the effects of ion-dust collisions and ion kinematic viscosities on the linear ion-acoustic 
instability as well as the nonlinear propagation of small amplitude solitary waves and shocks (SWS) 
in a negative ion plasma with immobile charged dusts. The existence of two linear ion modes, 
namely the 'fast' and 'slow' waves is shown, and their properties are analyzed in the collisional 
negative ion plasma. Using the standard reductive perturbation technique, we derive a modified 
Korteweg-de Vries-Burger (KdVB) equation which describes the evolution of small amplitude SWS. 
The profiles of the latter are numerically examined with parameters relevant for laboratory and space 
plasmas where charged dusts may be positively or negatively charged, ft is found that negative ion 
plasmas containing positively charged dusts support the propagation of SWS with negative potential. 
However, the perturbations with both positive and negative potentials may exist when dusts are 
negatively charged. The results may be useful for the excitation of SWS in laboratory negative 
ion plasmas as well as for observation in space plasmas where charged dusts may be positively or 
negatively charged. 
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I. INTRODUCTION 

The nonlinear propagation of solitary waves and shocks 
(SWS) in dusty plasmas have been widely studied for 
understanding the electrostatic disturbances in space 
plasma environments P, Q as well as in laboratory 
plasma devices 0, |j| ■ Many researchers have pointed out 
that charged dust grains can drastically modify the ex- 
isting response of electrostatic wave spectra in plasmas 
depending upon whether the charged dusts are consid- 
ered to be static or mobile [3l-fl2j . 

The existence of dust ion-acoustic (DIA) solitons was 
first predicted by Shukla and Silin @, and was later ob- 
served experimentally by Barkan et al The phase 
velocity of these DIA waves increases when the den- 
sity of electrons decreases. The nonlinear evolution of 
such small amplitude waves is described by the Korteweg 
de-Vries (KdV) equation, which was first derived by 
Washimi and Tanuiti [l3j in a two-component plasma. 
It was also reported that the KdV equation can have 
both compressive and rarefactive solitary wave solutions 
when sufficient amount of negative ions are present in 
multi-component plasmas [l4| . Furthermore, it was also 
observed that when the negative ion density exceeds a 
critical value, only negative or rarefactive solitons can be 
excited [HI]. 

One of the nonlinear phenomena of DIA waves is the 
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generation of shocks, which have already been observed 
experimentally (T(| . The shock waves are described with 
an additional Burger term due to dissipation in the KdV 
equation, known as the KdV-Burger (KdVB) equation 
Il6j . The nonlinear properties of dust ion-acoustic 
shocks and holes have been studied [l7| with the con- 
sideration of ion kinematic viscosity in a dusty plasma. 
Furthermore, the nonlinear properties of small amplitude 
ion-acoustic shocks in a collisional dusty plasma has been 
investigated by Ghosh et al [Isj . In their analysis, they 
have ignored the frictional force due to ion-dust collision 
and considered the effects of ion kinematic visc osity in 
dusty plasmas with positive ions. Mamun et al [111 ] in- 
vestigated the existence of dust electron-acoustic shocks 
in a negative ion plasma with dust charge fluctuation. 
They reported the formation of DIA shocks with negative 
potentials in the plasma. In a recent work of Adhikary 
fl2| , the effect of kinematic viscosity has been considered 
to study the propagation of DIA shocks in multi-ion plas- 
mas, and it has been stressed that the viscosity in dusty 
plasmas plays a key role in the formation of DIA shocks. 
Thus, both the viscosity and collision of ions with dusts 
may play crucial roles in dissipation on the propagation 
of SWS in dusty multi-ion plasmas. 

Since the charging rates of dusts by the positive and 
negative ions are nearly balanced, the existence and the 
formation of dusty plasmas with both positive and nega- 
tive ions may be possible under some circumstances. Re- 
cently, Kim et al li| have investigated that dust particles 
injected in laboratory negative ion plasmas can become 
positively charged when the number density of negative 
ions greatly exceeds (> 500) that of the electrons. In 
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space environments, the possible role of negative ions has 
been discussed, and it has been found that dusts can be 
positively charged if there is a sufficient number density 
of heavy negative ions (with mass > 300 amu) [2(| ■ 

In this work, we investigate the ion-acoustic instabil- 
ity as well as the nonlinear propagation of DIA SWS 
in an unmagnetized collisional dusty plasma containing 
both positive and negative ions. The effects of ion ther- 
mal pressures, ion-dust collisions as well as the kinematic 
viscosity of the ion fluids are taken into account to main- 
tain the equilibrium of ions. The present paper gener- 
alizes and extends the previous work of Adhikary [l2[ 
to include collisional effects of ions with dusts, different 
mass and thermal pressures of ions as well as different 
expression for ion kinematic viscosities. We show that 
when the negative ion concentration is much larger than 
the density of electrons or when plasma is the admixture 
of positively charged (static) dusts, SWS with negative 
potential exist. However, when the concentration of neg- 
ative ions is smaller than a certain value or plasma con- 
tains negatively charged dusts, both the compressive and 
rarcfactive SWS may exist. 



II. BASIC EQUATIONS 

We consider the one-dimensional propagation of DIA 
waves in an unmagnetized collisional dusty plasma, which 
consists of singly charged adiabatic positive and nega- 
tive ions, Boltzmann distributed electrons and immobile 
charged dusts. Since the dusts are too heavy to move 
on the time scale of the ion-acoustic waves, we do not 
consider the dynamics of charged dusts. The latter can, 
however, affect the collision rates with ions as well as the 
wave dispersion and nonlinearity. We have neglected the 
electron inertia, since the electron thermal speed is much 
larger than that of ions. We also assume that the nega- 
tive ions are heavier than the positive ions. The immobile 
dust particles carry some charges so as to maintain the 
overall charge neutrality condition given by 



n e o + n n0 = n p0 ± z d n d o, 



(1) 



where rijo is the unperturbed number density of charged 
species j (j = e, p, n, d, respectively, stand for electrons, 
positive ions, negative ions and static dusts), z d (> 0) 
is the dust charge state. The upper (lower) sign in Eq. 
([1]) corresponds to positively (negatively) charged dusts. 
The condition ([TJ can also be written as 



fl e = (M + (fl d - 1, 



(2) 



where fi e = n e0 /n n0 , fi { = n p0 /n n0 , fi d = z d n d0 /n n0 are 
the density ratios and £ = 1 (— 1) correspond to positively 
(negatively) charged dusts. The basic equations in one 
space dimension are 



drij d . . 



(3) 



dt + Ujd 



qj d(f> 3k B T 3 dn 2 ^ d 2 v 
1 Vj 



rrij dx 



2mjn/j dx 



dx 2 



\ (4) 



= 4-7re (n e - n p + n n - (z d n d0 ) , (5) 



and the Boltzmann distribution for electrons 



n e = n e0 exp (e4>/k B T e ) 



(6) 



where d/ dt = d t + Vjd x is the convective derivative. The 
physical quantities rij , Vj and rrij respectively denote the 
number density, velocity and mass of j-species particles. 
Furthermore, q p = e and q n = — e, where e is the elemen- 
tary charge. Also, (j> is the electrostatic potential, k B is 
the Boltzmann constant and Tj is the particle's thermo- 
dynamic temperature. In Eq. ((4]), Vj d (~ <Jj d n d QVtj) de- 
notes the collision rate in which <jj d is the collision cross 
section of j-species ions with dusts and Vtj = y/k B Tj /rrij 
is the thermal velocity of j-specics ions. Moreover, rfj 
is the kinematic viscosity of j-species ions which arises 
mainly due to the ion-dust collisions. We here mention 
that the exact expression for rfj in terms of Vj d or n d Q 
has not yet fully investigated so far. So, we will con- 
sider arbitrary values of it in order to see the damping of 
DIA waves numerically [2l| . We have used the adiabatic 
equation of state in Eq. (jU), namely Pj/Pjo — (nj/rijo) 3 
with Pjo 



ho 



k B T, 



lj for each of the ion species. The 
adiabatic index 7 = 3 [= (2 + D)/D, D being the num- 
ber of degrees of freedom] is due to the one-dimensional 
geometry of the system. Furthermore, in the ion con- 
tinuity equations the electron-ion recombination effects, 
being smaller at low-pressure limit (~ 4-5 x 10~ 4 Torr), 
have been neglected [16j |. However, the contribution from 
the frictional force Vj d Vj may not be neglected compared 
to that from the viscous term cc n, . It has been found 
that in laboratory experiments [2l|, the nondimensional 
viscosity parameter fjj = rfj/\ 2 D iu P n may be about 10 
times of the collisional frequency Vj d = i>j d /uj pn , where 
<jj pn — \J 47rn„oe 2 /m„ is the negative ion plasma fre- 
quency and Xd = \J k B T e / Ann n Qe 2 is the Debye length. 
Also, fjj and Vj d are almost linearly proportional to the 
dust number density n d Q. 

Next, we normalize the physical quantities according 



to (j) — > e(f>/k B T e 



rij/rijo, 



Vj/c s , where 



c s = ^/k B T e /m n is the ion-acoustic speed. The space 
and time variables are normalized by the Debye length 
Xu and the inverse of the negative ion plasma frequency 
u! pn respectively. Thus, from Eqs. ©-© we obtain the 
following normalized set of equations 



drij d . . 

-dT + dx M = > 



dt 



+ v id v 



Pi T 



dx 



3 dr% 
2 aj 8x 



- d 2y j 

Vj ~dx T 



(7) 
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d 2 4> 

dx 2 



(9) 



where the upper (lower) sign in =p on the right-hand side 
of Eq. © corresponds to positive (negative) ions, <7j = 
Tj/T e , pj = m n /mj with (3 p = (3, say. 



III. DISPERSION RELATION: ION-ACOUSTIC 
INSTABILITY 

In order to identify the different wave eigenmodes and 
to study their stability/ instability under the dissipation 
due to ion-dust collisions and ion kinematic viscosities, 
we perform a linear analysis in which the perturbations of 
the physical quantities vary as / ~ exp(ikx — iut), where 
ui (k) is the wave frequency (number). Thus, Fourier ana- 
lyzing the Eqs. 0-® we obtain the following dispersion 
relation 



D(u,k) = l+ ^ = < 



(10) 



where \ e = fi e /k 2 and 

X P = HiP [3/3a p k 2 - lo (w + w pd + ifjpk 2 )] , (11) 



Xn = \3a n k 2 - to 



ifj n k 2 )] 



(12) 



We note that the damping terms i(vjd + fjjk ) for j = p 
and n arise due to the ion-dust collisions and the viscos- 
ity effects, and thereby cause the DIA instability in the 
collisional plasma. When k 2 <C 1, fjjk 2 term is negligible 
compared to i)jd- Thus, for long wavelength modes, the 
damping effect is mainly due to the ion-dust collisions. 
In particular, in absence of these dissipative effects, the 
dispersion relation Eq. (fl"0|) reduces to [lj| 



Mi/3 



1 



lo 2 — 'Ak 2 j3(jp 



3£rer rl 



(13) 



The first and second terms on the right-hand side of Eq. 
(|13[) are, respectively, the contributions from the positive 
and negative ion species, whereas the term oc fx e is from 
the electron species in plasmas and the constant term 1 
is from the effect of dispersion due to charge separation 
(deviation from quasineutrality) of the species. If the 
phase speed of the wave satisfies u>/k S> c s > vtj , then 
from Eq. (fl"3]) the fast ion-acoustic wave mode can be 
recovered as 



(1 + M^)fc 2 



(14) 



In absence of negative ions in plasmas one obtains from 
Eq. (fT3")) the following disperson relation Q . 



2,(3a p k 2 



^k 2 



k 2 ' 



(15) 



From Eqs. (|T4|) and (jT5|) we find that for a fixed /3, /Xj 
and/or a p , the frequency of the ion wave modes increases 
with the wave number. However, as the wave number in- 
creases beyond a certain value, e.g., k > 0.2, the fast 
mode in the limit may approach more or less a constant 
value. Also, as the dust number density increases, i.e., fi e 
increases (decreases) for positively (negatively) charged 
dusts, the phase velocity of the fast mode decreases (in- 
creases). Physically, this decrease (increase) of the phase 
velocity is attributed to the addition (removal) of elec- 
trons by the dust grains in order to maintain the charge 
neutrality condition (fTJ) of the plasma. 

We numerically analyze the dispersion relation Eq. 
(|10p by considering parameters that may be represen- 
tative of laboratory negative ion plasmas with positively 
and negatively charged dusts. As in Ref. [lU, we con- 
sider a plasma in which positive ions are singly ionized 
K + and the heavy negative ions are SFq, such that 
/3 = m n /m p = 146/39 « 3.74. The temperature and 
density are considered as T e ~ T p ~ 0.2 eV, Th ~ T e /8 
with n„o ~ 2 x 10 9 cm~ 3 , n e o = n„o/700, n p o = 500n e o 
for positively charged dusts (and n e o = n„o/100, n p o = 
150n e0 for negatively charged dusts) to a surface poten- 
tial (j> s ~ 0.1 V and grain radius R = 5^m, so that 
Zd <~ i?0 s /e ~ 350. In this case, rido ~ 1.7 x 10 6 cm -3 
for which the rate of ion collisions with dusts would be 
higher than that with neutrals. We identify two ion wave 
modes, namely the 'fast' and 'slow' modes as exhibited 
in Figs. 1-3. 

Figure 1 shows the plots of the real wave frequency (up- 
per panel) and the growth rate (lower panel) of the DIA 
instability versus the wave number k obtained by numer- 
ically solving Eq. (fTOj) for positively charged dusts. The 
thick (black) line corresponds to the case when no dust 
is present in the plasma, whereas the thin (blue) line is 
for the case when dusts (positively charged) are present. 
In absence of ion thermal pressures, the real wave fre- 
quency and the growth rate assume almost a constant 
value for k > 0.2 as shown in Fig. 1 (red lines). This is 
evident from the analytic expression (|14|) . It is seen that 
the growth occurs, however, the frequency is reduced in 
presence of dusts. Thus, when dust is introduced into the 
plasma, the enhanced ion-dust collision frequency damps 
the instability rate. From Fig. 2, we observe that when 
the dusts are negatively charged, the real wave frequency 
and the growth rate become larger than that the case of 
no dust. This is in contrast to the case of positively 
charged dusts. The behaviors of these modes with or 
without charged dusts in plasmas suggest a possible di- 
agnostic for the roles of positively and negatively charged 
dusts as well as the effects of ion-dust collisions. 

Recently, it has been suggested by Rapp et al [2(j that 
the presence of sufficient amount (n„o > 40n e o) of heavy 
negative ions (with mass > 300 amu) may play important 
roles in the observation of positive dusts in the Earth's 
mesosphere between about 80 to 90 km. They also re- 
marked that the presence of negative ions corresponding 
to sub-nm size (0.3 — 0.4 nm) negatively charged dusts 




FIG. 1. (Color online) The Real (upper panel) and imaginary (lower panel) parts of uj versus k obtained as a numerical solution 
of the dispersion relation [Eq. QlQfr]. The thin (blue) and thick (black) lines correspond to the cases when positively charged 
dusts are present and when no dust is present in the plasma respectively. The line (red) which asymptotically approaches 
to more or less a constant value with k is corresponding to the case when ions are cold. The parameters are for laboratory 
negative ion plasmas with positively charged dusts and given by m n = 146m pro t, m p = 39m pro t (m pro t is the proton mass); 
T e ~ T p ~ 0.2 eV, T n ~ T e /8, n n0 ~ 2 x 10 9 cm" 3 , n e0 = n„o/700, n p0 = 500n e0 , n d0 ~ 1.7 x 10 6 cm" 3 , fj p = 0.2 and f] n = 0.1. 



may lead to the positive dusts with larger size (~ 2 nm) 
grains. In order to observe some features of the DIA 
modes in space plasmas, we consider plasma parame- 
ters that are representative of a dusty region at an al- 
titude of about 95 km where m n /m p = 300/28 w 10.7, 
T e ~ T p ~ T„ - 200 K, n n0 ~ 10 6 cm" 3 , n p0 ~ 2 x 10 4 
cm -3 and rido ~ 10 6 cm -3 . We assume that there is a 
population of sub-nm size negatively charged dusts in a 
dusty meteor trail (of radius R = 0.4 nm) region in the 
upper atmosphere. The features of the wave modes are 
shown in Fig. 3 for two different values of the density 
ratios: n n o/n e Q = 10 (thin line) and n n o/n e o = 50 (thick 
line). We find that the wave frequency and the growth 
rates are reduced by increasing the ratio. The qualitative 
behaviors of the other mode and the effects of charged 
dusts and the thermal pressure remain similar as Figs. 
1 and 2. In this case, we also note that the frequency 
and the growth rates increase faster with k having larger 
magnitudes compared to those in the case of laboratory 
plasmas described above. If such instabilities can occur, 
VHF/UHF radar scattering from DIA waves may be pos- 
sible diagnostic for the presence of sufficient amount of 
negative dusts in such meteor trail regions fl9l . l22j . 

It has been observed in experiments that slow wave 



modes do not favor the formation of solitons [15J, |23| . 
So, in the nonlinear regime, the fast mode may propa- 
gate as DIA solitary waves due to nice balance of the 
dispersive and nonlinear effects. However, as the dissipa- 
tion effects due to ion-dust collisions and ion kinematic 
viscosities, are entered into the dynamics, the system's 
evolution shows shock-like perturbations. These features 
will be studied in the next section. 



IV. DERIVATION OF THE EVOLUTION 
EQUATION 

We consider the nonlinear propagation of small but 
finite amplitude DIA waves in a collisional dusty negative 
ion plasma. In order to observe the evolution of the waves 
in a different frame of reference which is moving with a 
speed vq, we stretch the coordinates as £ = e 1 ' 2 (x — vot) 
and t = e 3 / 2 £, where e is a small parameter measuring the 
weakness of perturbations. We also assume that Vjd = 
e 3 / 2 i>jo and rjj = e 1 ^ 2 rjjo, where Vj$ and 77^0 are of the 
order of unity or less. Such a consideration of smallness 
of T)j and Vjd can be found in the literature [l8j . and 
is valid in many experimental situations (see e.g., Ref. 
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FIG. 2. (Color online) The Real (upper panel) and imaginary (lower panel) parts of uj versus k obtained as a numerical solution 
of the dispersion relation [Eq. (fTtJj) ]. The thin (blue) and thick (black) lines correspond to the cases when negatively charged 
dusts are present and when no dust is present in the plasma respectively. The line (red) which asymptotically approaches to 
more or less a constant value with k is corresponding to the case when ions are cold. The parameters are for laboratory negative 
ion plasmas with negatively charged dusts and given by n„o ~2x 10 9 cm" 3 , n e o = n n o/100, n p o — 150n e o and ndo ~ 2.8 x 10 6 
cm -3 . Other parameters are the same as in Fig. 1. 



|21|). The dynamical variables are expanded as 



1 + enf ] + e 2 nf ) + 



v j =evf ) +<?vf> + 
6 = e6^ +e 2 </>( 2) + • 



(16) 



(17) 



We then substitute these expansions and the stretched 
coordinates into Eqs. ©-©j and equate different powers 
of e. In the lowest order of e (i.e., e 3 / 2 ) we obtain the 
following relations for the first order perturbations 



(i) 

n) = a j 



(18) 



together with the expression for the wave speed in the 
moving frame of reference as 



1 

2jU e 



S ± \S 2 - 12/3^e [cr P + cr„(/i, + 3^ e o>)] 



(19) 

where 5=1 + /ij/3 + 3(cr„ + j3a p )iJ, e and ctj = ±/3j/ {vq — 
3j3j(jj). Here ± denotes the quantities corresponding to 
positive (J = p) and negative (j = n) ions respectively. 
This ± sign appears due to the consideration of differ- 
ent mass and temperatures of positive and negative ions. 



However, in a limited situation where m p = m n and 
T p =T n , we can recover the expression (here one should 
consider the plus sign before the square root in Vq) as in 



Ref. [12(. However, in the next section, we will see that 



only the minus sign (of ±) favors the formation of SWS in 
laboratory and space plasmas, since for the plus sign, the 
nonlinear term becomes much larger than the dispersive 
term in the evolution equation to be shown shortly. 

Proceeding to the next order of e (i.e., e 5 / 2 ) we ob- 
tain the following set of equations for the second order 
perturbed quantities 



a ( 2 

UTl j ' 

- Vo —^7 h Oi j 



da 



~d7 



2 h^y 

■CtjVo- 



dv 



(2) 



0, (20) 



0_ 



+ PjoJ a-jvocp^ + i« 2 (v 2 + 2>Pj(jj) 




(21) 
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FIG. 3. (Color online) The Real (upper panel) and imaginary (lower panel) parts of u versus k obtained as a numerical 
solution of the dispersion relation [Eq. (|10[> ] . The thin (red) and thick (black) lines correspond to the two values of the density 
ratio: n n o/n e o — 10 and 50 respectively. The parameters are for space plasmas with negatively charged dusts and given by 
m n — 300m prot , m p = 28m pr ot (m pro t is the proton mass); T e ~ T p ~ T n ~ 200 K, n„o ~ 10 6 cm" 3 , n p o = 60n e o and n d0 ~ 10 6 
cm" 3 . 




(22) 



where ± sign in Eq. (|21|) corresponds to positive (j = p) 
and negative (J = n) ions respectively. 

Eliminating the second order quantities from Eqs. 
(f20)) - ([2"2")) . we obtain, after few steps, the following evolu- 
tion equation of the Korteweg de-Vries-Burger (KdVB) 
type. 



<9$ , <9$ <9 3 $ 



d 2 ® 



(23) 



where $ = 4>^ x \ The coefficients of nonlinearity and dis- 
persion as well as of dissipation due to the ion kinematic 
viscosities and the ion-dust collisions are, respectively, 
given by 



.4 



3a^i(vl + /3a p ) + 3(3al(v 2 + a n ) + (3(a n - ^a p ) 



D 



(24) 
(25) 



2(^1 + f3a 2 n ) 



2{ma 2 p + Pa 2 n ) 



(26) 



In particular, for 



= mp, T n = T„ and r/ n = r) p , we 
recover the expressions as in Ref. [12(. When the dissi- 
pation terms proportional to r\ and v are ignored, i.e., in 
absence of frictional force and viscous stress, the resul- 
tant equation is the KdV equation discussed in Ref. . 
Thus, the Burger term in Eq. (|23|) gives rise the gen- 
eration of shocks and the frictional force term provides 
damping of the wave due to ion-dust collision. The latter 
increases with the increase of dust number density. 



V. RESULTS AND DISCUSSION 

Typically, the existence of SWS with positive or nega- 
tive potential depends on the sign of the nonlinear coeffi- 
cient A. So, we numerically examine the behaviors of the 
nonlinear and the dispersion coefficients A and B with 
the variations of the positive to negative ion density ra- 
tio /!.;. The parameters are considered for both positively 
(Fig. 4) and negatively (Fig. 5) charged dusts as in Fig. 
1. The solid and dashed lines correspond to T n = T e /8, 
T p - T e and T n = T e /10, T p - T e /2 respectively. We 
have considered, respectively, the negative and positive 
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signs in the expression for Vq [Eq. (p~9|) ] for the subplots 
[(a), (c)] and [(b), (d)] of Figs. 4 and 5 to compare the 
corresponding values of A and B. From Figs. 4(b) and 
(d) we find that for multi-ion plasmas with positively 
charged dusts, the numerical values of the dispersion co- 
efficient B are much larger than those of the nonlinear co- 
efficient A. However, for plasmas with negatively charged 
dusts, the coefficient A can become larger than B with an 
increasing value of [See subplots 5(a) and (c)]. So, in 
both these cases, the balance between the nonlinear and 
the dispersive effects may not occur for the formation 
of KdV solitons in collisional dusty negative ion plas- 
mas. Furthermore, when the dissipation overwhelms the 
wave dispersion and it is in nice balance with the nonlin- 
earity arising from the nonlinear mode coupling of finite 
amplitude waves, the monotonic shocks may be gener- 
ated. It turns out that the values of A and B in the 
subplots [(a), (c)] of Fig. 4 may favor the formation of 
SWS with negative potentials (Since A maintains only 
the negative sign), whereas those in the subplots [(b), 
(d)] of Fig. 5 may favor the formation of SWS with both 
positive and negative potentials (Since A can have both 
the positive and negative signs). Thus, we may conclude 
that plasmas with positively charged dusts may support 
the propagation of SWS with negative potentials, how- 
ever, SWS with both positive and negative potential may 
exist when dusts are negatively charged. In the former 
case, the number density of negative ions greatly exceeds 
that of electrons and positive ions, whereas, positive ion 
density is to be larger than the density of negative ions in 
the latter case. In our numerical solution of Eq. (j2"3")l be- 
low, we present the characteristic features of those SWS 
which have negative potential (as the qualitative features 
of SWS with positive potential will be similar). 

Next, we numerically solve Eq. ([23]) by Runge-Kutta 
scheme with an initial condition of the form $(£) = 
— 0.03scch 2 (£/15). The development of this waveform at 
different times and for different parameters are shown 
in Figs. 6-8. The parameter values are considered as 
in Figs. 1-3 for both positively and negatively charged 
dusts in electronegative plasmas. Figure 6 shows the case 
of KdV soliton, i.e., in absence of any dissipation for dif- 
ferent times: (a) r = 0, (b) r = 200, (c) r = 300 and 
(d) r = 500. We see that the leading part of the initial 
pulse steepens due to positive nonlinearity. As the time 
progresses, the pulse separates into solitons and a residue 
due to the wave dispersion. It is clear from Fig. 6(d) that 
once the solitons are formed and separated, they propa- 
gate without changing their shape due to the nice balance 
of the nonlinear and dispersion terms. This subplot also 
shows that (see the dashed line) when charged dusts are 
introduced into the plasma, the magnitude of the soliton 
height is reduced, but the width is increased. However, 
as the number density of negative ion increases (see the 
dotted line), both the amplitude and width of the soliton 
decrease. In each of these cases the soliton is seen to be 
up-shifted. 

Figure 7 shows the profiles of the solution of KdVB 



equation for plasmas with positively charged dusts. In 
this case, the values of rjjo and Vjo are considered arbi- 
trarily in order to observe the damping of solitary waves 
[2l] |. From Fig. 7(b), it is clear that as the dissipation 
due to ion kinematic viscosity enters into the dynam- 
ics, the solitary pattern breaks up and a shock profile is 
formed. Also, as the term oc rjj increases, the height of 
the shock profile decreases. Physically, the increase of 
rjj corresponds to the increase of the dust number den- 
sity, which, in turn, increases the ion-dust collision fre- 
quency. While compared with Fig. 7 (a), Figures 7(c) 
and (d) show the wave damping due to the presence of 
ion-dust collisions. The latter are shown to decrease the 
wave amplitude significantly as in Fig. 7(d) when some 
larger values of them are considered. From Figure 7(c) 
we find that the presence of charged dusts as well as the 
increase of the negative ion density decrease the ampli- 
tude of shocks. This is expected as the introduction of 
charged dusts causes dissipation via the collisional term 
and the viscosity term. Also, as the negative ion density 
increases, more electrons will be removed by the dust 
grains in order to maintain the charge neutrality. This 
also results to the increase of the positively charged dust 
density and hence may cause the reduction (in magni- 
tude) of the wave amplitude. 

Figure 8 shows the developments of solitary waves into 
shocks for plasmas with negatively charged dusts. In ab- 
sence of the dissipative effects [Fig. 8(a)], a different wave 
pattern with a train of solitons is seen to form due to the 
dominant role of wave dispersion over the nonlinearity. 
As the dissipation due to viscosity effects starts playing 
a role, the number of wave trains reduces and the wave 
steepening begins to occur at a value of rjj < 1 [See Fig. 
8(b)]. Further increasing the value of rjj (rjj ~ 1) results 
into the fact that the train of solitons disappear leaving 
only one wave front with a lower value of |$| [See Fig. 
8(c)]. From Fig. 8(d) it is seen that due to the effects of 
ion-dust collisions, the wave gets damped [Similar to the 
case of positively charged dusts (Fig. 7)] and the wave 
amplitude reduces in its magnitude. We note that the 
damping seen in this case is qualitatively different from 
the case of positively charged dusts (Fig. 7). From Fig. 
8(d) we also find that in contrast to Fig. 7(c) where the 
negative ion density decreases the wave amplitude, as the 
number density of positive ions increases, the magnitude 
of $ increases. Here the increase of the positive ion den- 
sity causes the addition of more electrons into the dust 
grain surface, which, in turn, increases the dust charge 
state in order to maintain the charge neutrality. As a re- 
sult, the dissipative effects due to dust density enhance- 
ment become more pronounced and the wave steepening 
occurs with an increased |$|. 

Another interesting feature may be the formation of 
double-layer solutions of Eq. (|2"5j) . Such solutions con- 
sist of two layers with positive and negative space charge 
and support localized electric fields. Also, double-layers 
(DLs) can play important roles in providing a mechanism 
for supporting electric fields in collisionlcss plasmas with 
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FIG. 4. (Color online) The plots of the nonlinear coefficient A (upper panels) and the dispersion coefficient B (lower panels) 
versus the positive to negative ion density ratio fii are shown for plasmas with positively charged dusts. The other parameters 
are as in Fig. 1. The solid and dashed lines correspond to T n = T e /8, T p ~ T e and T n = T e /10, T v ~ T e /2 respectively. In the 
subplots [(a), (c)] and [(b), (d)], the values of Vq [Eq. <|19p ] corresponding respectively to the negative and positive signs have 
been considered to compare the results as well as to show that the numerical values of A and B in the subplots [(a), (c)] may 
favor (while the others may not) the formation of solitary waves and shocks with negative potentials. 



nearly zero resistivity. Furthermore, DLs have relevance 
to charge particle acceleration in cosmic plasmas and in 
plasma thrusters. However, there are several mechanisms 
by which DLs can be formed (See, e.g., Ref. [24| for some 
recent review works on DLs). Nevertheless, such DLs are 
more difficult to generate in a laboratory and require a 
fine tuning of the plasma parameters. However, the detail 
discussion on the formation and the properties of DLs is 
limited to the present study as we have mainly focused on 
the formation and the characteristics features of solitons, 
and shocks so generated due to dissipative effects. 

In what follows, Eq. (|23|) can have a travelling wave 
solution, in absence of the collisional effects (y = 0), 
given by po| 



$(£,t) = Ci 



^ = -— z 



12-q 2 



25AB(1 
\ 5B 



125AB 2 



(27) 



where C\ and Ci are arbitrary constants. The profiles 
of Eq. (f27|) for different plasma parameters are shown 
in Fig. 9. Figures 9(a)-(c) exhibit some shock profiles 
corresponding to plasma parameters relevant for nega- 
tively charged dusts, whereas Fig. 9(d) is for negative 



ion plasmas with positively charged dusts. We find that 
the DLs with both the compressive and rarefactive po- 
tentials may exist [Figs. 9(a) and (b)] for a certain range 
of values of r]j , beyond which DLs with only negative 
potential (rarefactive) may be formed [Fig. 9(c)]. From 
Figs. 9(a) and (b), it is also seen that for a fixed negative 
ion density as the density of positive ions increases, the 
value of | | increases by the same physical reason as in 
Fig. 8(d). The opposite trend may be seen to occur by 
increasing the negative ion density (in case of positively 
charged dusts) or by decreasing the values of r]jo . The 
effects of the latter are shown in Fig. 9(c) exhibiting DLs 
with only negative potential. It turns out that the wave 
front with two layers may not be formed by gradually de- 
creasing the values of r]j or by reducing the charged dust 
state or density in negative ion plasmas with negatively 
charged dusts, rather we can see the shock wave fronts 
with only negative potentials. Figure 9(d) also shows 
that if the dusts are positively charged or electrons are 
almost absorbed by the dust grains in pair-ion plasmas, 
the formation of DLs may not be possible. However, a 
shock profile with negative potential may exist with a 
significant drop of the wave amplitude |<I>| compared to 
the case of negatively charged dusts. 
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FIG. 5. (Color online) The plots of the nonlinear coefficient A (upper panels) and the dispersion coefficient B (lower panels) 
versus the positive to negative ion density ratio fj,i are shown for plasmas with negatively charged dusts. The other parameters 
are as in Fig. 1. The solid and dashed lines correspond to T n = T e /8, T p ~ T e and T n = T e /10, T p ~ T e /2 respectively. In the 
subplots [(a), (c)] and [(b), (d)] the values of Vq [Eq. (|I9|l ] corresponding respectively to the negative and positive signs have 
been considered to compare as well as to show that the numerical values of A and B in the subplots [(b), (d)] may favor (while 
the others may not) the formation of solitary waves and shocks with both positive and negative potentials. 



To summarize, we have investigated the propagation 
characteristics of DIA waves in a collisional negative ion 
plasma with immobile charged dusts. The latter may be 
positively charged when the number density of negative 
ions exceeds that of positive ions and is much larger than 
that of electrons. In the linear regime, we have studied 
numerically the ion-acoustic instability with plasma pa- 
rameters relevant for both laboratory [22[ and space 
plasmas [20[. We find that the two modes, namely 'fast' 
and 'slow' waves exist in dusty negative ion plasmas. Due 
to the higher values of the phase velocity of the fast ion 
wave than the ion thermal speeds, the Landau damping 
effect on dissipation is negligible, however, the damping 
of the solitary waves is mainly caused by the ion-dust 
collisions and the ion kinematic viscosities. We find that 
the slow modes propagate without any instability and 
with frequency below the negative ion plasma frequency. 
However, for long wavelength fast modes, the damping 
effect is mainly due to the ion-dust collisions. For lab- 
oratory negative ion plasmas when the electron flow is 
sufficient to drive the ion-acoustic instability, the inclu- 
sion of micron seized dust grains may damp the instabil- 
ity due to their collision with ions. The frequency and 
the growth rates are reduced when dusts are positively 



charged. However, the opposite trend occurs when dusts 
are negatively charged. Thus, the phase velocity of the 
waves may decrease or increase with the dust density ac- 
cording to when the dusts are positively or negatively 
charged. These results suggest a possible dilagnostic 
for the presence of positively or negative charged dusts 
as well as for the effects of ion-dust collisions in dusty 
negative ion plasmas [l9|, [22[. For the space situation, 
though the behaviors remain similar, however, the wave 
frequency and the instability rate increase faster with the 
wave number than the case of laboratory plasmas. 

We have also studied the nonlinear propagation of 
small amplitude fast modes as DIA SWS in dusty neg- 
ative ion plasmas. We show that the evolution of such 
waves can be described by a modified KdVB equation. 
The latter generalizes and modifies the previous investi- 
gation in negative ion plasmas [l2| . The KdVB equation 
is numerically solved to show that the perturbations with 
negative potential may propagate as SWS in plasmas 
with positively charged dusts, whereas SWS with both 
positive and negative potential may exist when dusts are 
negatively charged. The profiles of these SWS with only 
negative potential are shown graphically (as those with 
positive potential are similar), and analyzed with pa- 
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FIG. 6. (Color online) Development of an initial pulse of the form $(£) = — 0.03sech 2 (£/15) into solitary waves at different 
times for the parameters as in Fig. 1, and when }/jd = fjjd = 0. In the subplot (d), the dashed (red) and dotted (black) lines 
correspond to the cases of no dust and when negative ion density is different, i.e., n„o = 1000n e o in the plasma respectively. 



rameters relevant for laboratory and space plasmas. We 
find that, in contrast to the effects of negatively charged 
dusts, the presence of positively charged dusts reduces 
the wave amplitude, but enhances the width of solitary 
waves. The theoretical results may be useful for the ob- 
servation of dust ion-acoustic waves in space plasmas, 
e.g., a dusty meteor trail region in the upper atmosphere 
as well as the experimental verification of the excitation 
of ion-acoustic instability and the nonlinear propagation 
of ion-acoustic solitary and shock waves in dusty multi- 
ion plasmas. 

There may be some open issues, e.g., the theoretical 
deduction of ion kinematic viscosity and the effects of 
ion drag forces due to positive and negative ions on the 
dust particles could be problem of interest but beyond the 



scope of thel present investigation. Such forces, which 
are opposite in direction, may also be comparable in mag- 
nitude when the positive and negative ion densities arc 
comparable. Since the outward ion drag force is known 
to be responsible for the formation of a void under mi- 
crogravity conditions, the presence of both the positive 
and negative ions could forbid the formation of voids. 
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FIG. 8. (Color online) Development of an initial pulse of the form $(£) = — 0.03sech 2 (£/15) into shocks at r = 500 for different 
values of Vjo and r/jo as in the figure. The parameters are for plasmas with negatively charged dusts where n e o = n„o/10, 
n p o — 20n e o. Others parameters are as in Fig. 1. In the subplot (d), the dashed (black) line corresponds to different value of 
the positive ion density, i.e., n p o = 40n e o with the same other parameters as the solid line. 



13 




FIG. 9. (Color online) Analytic solution of Eq. (|27[) for different plasmas with negatively [subplots (a)-(c)] and positively 
charged dusts [subplot (d)]: (a) n e0 = n„ /2, n p0 = 5n e0 , r] p0 = 1.2, rj n0 = 1.1; (b) n e0 = n„ /2, n p0 = 10n e o, VpO = 1-3, 
r] n0 = 1.2; (c) n e o = n n0 /2, n p0 = 5n e0 , 77 P o = 0.7, 77„ = 0.6; (d) n e0 = n n0 /700, n pa = 500n e o, r) pa = 1.2 and rj„ = 1-1. 



